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Using advanced lattice methods in Quantum Chrome-dynamics, three distinct scales are estab- 
lished in the excitation spectrum of the gluon field around a static quark-antiquark pair as the 
color source separation R is varied. On the shortest length scale, the excitations are consistent with 
states created by local gluon field operators arising from a multipole operator product expansion. 
An intermediate crossover region below 2 fm is identified with a dramatic rearrangement of the 
level orderings. On the largest length scale of 2-3 fm, the spectrum agrees with that expected for 
string-like excitations. The energies nearly reproduce asymptotic tt/R string gaps, but exhibit a 
fine structure, providing important clues for developing an effective bosonic string description. 

PACS numbers: 11.15.Ha, 11.25.Pm, 12.38.Aw, 12.38.Gc 



There exists great interest and considerable effort to 
explain quark confinement in Quantum Chromodynamics 
from the string theory viewpoint. The ideas of 't Hooft, 
Polyakov, Wittcn, and others, and recent glueball spec- 
trum and string tension calculations in AdS theories are 
illustrative examples. In a somewhat complementary ap- 
proach, the search for a microscopic mechanism to ex- 
plain quark confinement in the QCD vacuum contin- 
ues with vigorous effort. Our view is that the relevant 
properties of the underlying effective QCD string theory, 
whether it emerges from strictly string theoretic ideas or 
from the microscopic theory of the vacuum, are coded in 
the excitation spectrum of the confining flux. To estab- 
lish the main features of this spectrum is the objective of 
this work. 

We report here, for the first time, our comprehen- 
sive determination of the rich low-lying energy spec- 
trum of the excited gluon field between a static quark- 
antiquark (qq) pair in the fundamental color represen- 
tation of SU(3) C . Only a few of these states have been 
previously studied and only for qq separations R up to 1 
fm 0, 0) Q ■ Due to our use of recent advances in lattice 
gauge theory technology, which include anisotropic lat- 
tices, an improved lattice gauge action, and large sets of 
creation operators, we have been able to determine the 
excitation spectrum for a large range of R values from 
0.1 to 3 fm, and for a greatly extended set of excited 
states never studied before. To investigate the onset of 
string-like behavior when the flux is much longer than 
its intrinsic width became a realistic goal Q . The effects 
of light qq pairs in the vacuum are not included in this 
study, with focus on the confining properties of the gluon 
field. Earlier results from less extensive simulations were 
reported in Ref. [5j. 



Three exact quantum numbers which are based on the 
symmetries of the problem determine the classification 
scheme of the gluon excitation spectrum in the presence 
of a static qq pair. We adopt the standard notation from 
the physics of diatomic molecules and use A to denote 
the magnitude of the eigenvalue of the projection J g ■ R 
of the total angular momentum 3 g of the gluon field 
onto the molecular axis with unit vector R. The capital 
Greek letters A, $, ... are used to indicate states 

with A = 0, 1, 2, 3, . . ., respectively. The combined oper- 
ations of charge conjugation and spatial inversion about 
the midpoint between the quark and the antiquark is also 
a symmetry and its eigenvalue is denoted by rjcp ■ States 
with r\cp = 1( — 1) are denoted by the subscripts g (u). 
There is an additional label for the S states; S states 
which are even (odd) under a reflection in a plane con- 
taining the molecular axis are denoted by a superscript 
+ (— ). Hence, the low-lying levels are labeled £+, E~ , 



A g , A„, and so on. For convenience, 



we use r to denote these labels in general. 

The gluon excitation energies Er(R) were extracted 
from Monte Carlo estimates of generalized large Wilson 
loops using extended sets of 15 to 25 operators, projected 
onto the various symmetry sectors, on three-dimensional 
time slices of the lattice. Monte Carlo estimates of the 
Wy (R, t) Wilson correlation matrices were obtained in a 
series of simulations performed on our dedicated UP2000 
Alpha cluster using an improved gauge-field action 
The couplings /3, input aspect ratios £, and lattice sizes 
for each simulation will be detailed in forthcoming pub- 
lications but are available now upon request. The use of 
anisotropic lattices, with the temporal lattice spacing at 
much smaller than the spatial spacing a s , was crucial for 
resolving the gluon excitation spectrum, particularly for 
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large R. The renormalized anisotropy was determined us- 
ing the torelon dispersion relation. To hasten the onset of 
asymptotic behavior in effective mass plots, iteratively- 
smeared spatial links were used in the generalized Wilson 
loops whose temporal segments were constructed from 
thermally-averaged links to reduce statistical noise. 

Restricted to the R = 0.2 — 3 fm range of a selected 
simulation, energy gaps AE above the ground state are 
compared to asymptotic string gaps for 18 excited states 
in Fig. □ The quantity AE/(Ntt/R) - 1 is plotted to 
show percentage deviations from the asymptotic string 
levels for string quantum number N = 1,2,3,4. The 
energy gaps, far below the null lines of the plots and 
strongly split at fixed N, differ from the simple string 
gaps for R < 2 fm and a fine structure remains visible 
beyond R = 2 fm. The three relevant scales of the spec- 
trum are discussed below. 

The errors shown in Fig. 1 are statistical only and were 
determined using standard bootstrap techniques. Au- 
tocorrelations were explicitly calculated to ensure that 
the statistical errors were not underestimated. Depen- 
dence of the error estimates on the choice of temporal fit 
ranges was minimized by using fits to sums of two ex- 
ponentials. The removal of contamination in our results 
from excited-states was checked using a variety of meth- 
ods, including Bayesian constrained fits to large sums of 
exponentials, variational re-diagonalizations, and corre- 
lation matrix fits to sums of coupled exponentials. The 
insignificance of systematic effects from finite volume was 
verified, and simulations carried out for several other lat- 
tice spacings indicate that finite spacing errors in the 
results shown in Fig. 1 are expected to be small. 
Short distance spectrum. For R <C 1 fm, the ob- 
served level ordering is consistent with the short-distance 
physics of gluon field excitations which are trapped 
around a dominantly color octet static qq pair. Since 
gluon field dynamics will only depend softly on R, a 
new approximate symmetry is expected in the R\[5 <C 1 
limit (the string tension a sets the scale of transverse di- 
mension across the flux). This symmetry was first noted 
within the context of the static bag picture of gluon exci- 
tations @ and used more extensively in Refs. pa. 15 1^. Hfl| . 
In particular, it was shown in Ref. 10] how the short- 
distance operator product expansion of gluon excitations 
around static color sources, in lowest order of the multi- 
pole expansion of the gluon field, can be applied to the 
spectrum yielding similar results to the bag picture |£| 
without phenomenological model assumptions. 

Gluon excitations will transform according to repre- 
sentations of 0(3) in the short-distance limit so that A 
in the set T will be replaced by the gluon field angular 
momentum J. A representation with fixed J contains 
degenerate states with different A values. For illustra- 
tion, gluon field operators of 0(3) multiplets are listed in 
Table [|] for ten states with their short distance angular 
momentum quantum number J for each multiplet. There 
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FIG. 1: Energy gaps AE above Ejj" are shown in string units 
for quantum numbers in continuum and lattice notation. The 
Nambu-Goto string is discussed in the text. 



are two triplets (£+ , n g ), (E u , II U ), and two J = 2 mul- 
tiplets (Sg ,n g , A g ), (E+,l4,A u ) in TableD II and A 
states occur with two opposite sign chiralities A within 
each 0(3) multiplet. Within the degenerate multiplets, 
one finds different string quantum numbers mixed to- 
gether. Fig. [3 illustrates the remarkable working of the 
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TABLE I: Operators to create excited gluon states for small 
qq separation R are listed. E and B denote the electric and 
magnetic operators, respectively. The covariant derivative D 
is defined in the adjoint representation 10]. 



gluon state 
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operator 
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R E, R (D x B) 


n 9 
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predicted short-distance degeneracies. Only the states 
A u and £+ show considerable soft breaking of the ap- 
proximate symmetry at the shortest R values. 
Crossover region. For 0.5 fm < R < 2 fm, a dramatic 
crossover of the energy levels toward a string-like spec- 
trum as R increases is observed. For example, the states 
£~ with N = 3 and £~ with N = 4 break violently away 
from their respective short-distance 0(3) degeneracies to 
approach the ordering expected from bosonic string the- 
ory near R ~ 2 fm. 

An interesting feature of the crossover region is the suc- 
cessful parametrization of the ground state energy by 
the empirical function Eq (i?) = a + <tR — , with the 
fitted constant c close to unity, once R exceeds 0.5 fm. 
The Casimir energy of a thin flux line was calculated in 
Refs. 0,0] , yielding c = 1, and this approximate agree- 
ment is often interpreted as evidence for string formation. 
While the spectrum, including the qualitative ordering 
of the energy levels, differs from the naive bosonic string 
gaps for R < 1 fm, a high precision calculation shows 
the rapid approach of c c g to the asymptotic Casimir 
value in the same R range |l3j . Although there is no in- 
consistency between the two different findings, a deeper 
understanding of this puzzling situation is warranted. 

We will return to this issue in a high precision study of 
the 3-dimensional Z(2) gauge model in a future publica- 
tion Q. This accurate study of c g(R) and the excita- 
tion spectrum of the Z(2) flux line for a wide range of R 
values between 0.3 fm and 10 fm will clearly demonstrate 
the early onset of c « 1 without a well-developed string 
spectrum. For now, Fig.|3|shows the lowest excitations in 
Z(2) for R = 0.7 fm, revealing a bag-like disorder profile 
surrounding the static qq pair in the vacuum ^ij. The 
two lowest energy levels are substantially dislocated from 
exact it / R string gaps and all other excitations form a 
continuous spectrum above the glueball threshold. Since 
the submission of this work, a new study of Z(2) at fi- 
nite temperature has appeared |l5| . reporting very early 
onset of string behavior in support of Ref. |13| . 
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FIG. 2: Short-distance degeneracies and crossover in the 
spectrum. The solid curves are only shown for visualization. 
The dashed line marks a lower bound for the onset of mixing 
effects with glueball states which requires careful interpreta- 
tion. 

String limit. For R > 2 fm, the energy levels exhibit, 
without exception, the ordering and approximate degen- 
eracies of string-like excitations. The levels nearly re- 
produce the asymptotic n/R gaps, but an intriguing fine 
structure remains. 

It has been anticipated that the interactions of mass- 
less excitations on long flux lines are described by a lo- 
cal derivative expansion of a massless vector field £ with 
two transverse components in four-dimensional space- 
time 0, 0]. Symmetries of the effective QCD string 
Lagrangian require a derivative expansion of the form 

(1) 

where the dots represent further terms with four or more 
derivatives in world sheet coordinates. The coefficient a 
has the dimension of a mass squared and can be identified 
with the string tension a. The other coefficients must be 
determined from the underlying microscopic theory. Ex- 
amples with calculable coefficients include the D — 3 Z(2) 
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FIG. 3: The disorder profile of the vacuum and two excita- 
tions are shown with arrows indicating the qq axis. 



defect line and the D = 4 Nielsen- Olesen vortex in a semi- 
classical expansion of field theory. Similar calculations 
in QCD would require a quantitative understanding of 
quark confinement. 

Asymptotic string gaps and the observed fine struc- 
ture of the spectrum provide the most direct tests of 
effective QCD string formation for large R. Low fre- 
quency excitations are dominated by the first term which 
describes the asymptotic spectrum of massless modes 
with exact tt/R string gaps. Higher dimensional oper- 
ators are expected to generate perturbative fine struc- 
ture corrections in higher powers of Matching 
the coefficients of Eq. (1) to the spectrum is a real- 
istic goal for future high precision simulations. Since 
it is reasonable to expect that the first few coefficients 
arise from the geometric properties of the effective string, 
it is useful to introduce the Nambu-Goto (NG) action 
cr-y/1 + <9 M £ ■ <9 M £ as the simplest geometric model. The 
NG spectrum with fixed end boundary conditions in D 



dimensions was first calculated in Ref. 16] with the result 
E N = aR(l - -§^tt + fjp-)^- Although there exists 
a quantization problem in angular momentum rotations 
around the qq axis at finite R values unless D = 26 0] , 
the problem disappears in the large R limit. It has been 
widely expected that predictions at large R from the NG 
model will be similar to effective bosonic string theory 
with a good description of the observed spectrum. De- 
viations from the NG string at large R in Fig. ^ ma y 
indicate rigidity and other missing physical properties of 
the flux line in the model. 

In this work, striking confirmation of string-like flux 
formation of the gluon field for large qq separations was 
presented. Our simulations revealed a first glimpse of a 
fine structure in the large-distance spectrum, offering the 
tantalizing possibility of understanding the nature of the 
QCD string in future high precision simulations. 

Throughout this extended project, we benefited from 
discussions with M. Peardon, G. Bali, R. Brower, 
K. Intriligator, P. Hascnfratz, M. Liischer, C. Michael, 
S. Necco, F. Niedermayer, R. Sommer, J. Soto, P. Weisz, 
and many others. This work was supported by the DOE, 
Grant No. DE-FG03-97ER40546, the NSF under Award 
PHY-0099450 and the European Community's Human 
Potential Programme under contract HPRN-CT-2000- 
00145, Hadrons/Lattice QCD. 
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